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METHODS AND APPARATUS FOR ENHANCING DIAGNOSIS OF 

MYOCARDIAL INFARCTIONS 

CROSS-REFERENCE TO RELATED APPLICATION 
[0001] This is a regular patent application claiming the benefit under 35 U.S.C. §1 19(e) from 
U.S. Provisional Patent Application No. 60/396,681 filed July 17, 2002, the full disclosure of 
which is incorporated herein by reference. 

[0002] Various embodiments of the present invention employ methods, apparatus, and/or 
systems related to similar methods, apparatus, and/or systems described in PCT Application No. 
WO 01/67954 Al, filed March 14, 2001, by a common inventor to the present application, the 
entire contents of which are hereby incorporated by reference. A copy of WO 01/67954 Al is 
attached to this application as Appendix A. 

BACKGROUND OF THE INVENTION 
[0003] The present invention relates generally to diagnosis of myocardial infarctions. More 
specifically, the invention relates to methods and apparatus for enhancing diagnosis and 
localization of acute myocardial infarctions using electrocardiogram data. 

[0004] Acute myocardial infarction (AMI), commonly known as heart attack, is the single 
most common cause of death in the US. This disease entity is typically caused by an acute 
thrombotic event which leads to occlusion of a major coronary artery. When a coronary artery is 
occluded (i.e., blocked), the portion of the heart tissue which that artery supplies with blood no 
longer receives sufficient blood to provide oxygen to the tissue. If not quickly treated, this lack 
of blood supply ("ischemia") results in permanent damage and cell death in the myocardial 
tissue. Since myocardial damage due to cardiac ischemia becomes irreversible after a certain 
time, it is imperative that therapy to restore blood flow to the heart be started as quickly as 
possible. Thus, rapid diagnosis of AMI is also imperative, so that appropriate treatment can be 
started as quickly as possible. 

[0005] Typical therapies for AMI include administering thrombolytic agents intravenously or 
percutaneous transluminal coronary angioplasty (PTCA) with or without stent placement. In 
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either type of therapy, it would be advantageous to have information regarding the location of the 
AMI on the heart and the identity and/or location of the coronary artery that is blocked. This 
would allow the medical professional administering the therapy to focus quickly on the affected 
area for treatment. Information about where on the heart the AMI has occurred and information 
about the size of the AMI would help a medical professional to determine which artery is 
occluded and where, along the length of the artery, the occlusion is located. 

[0006] The most commonly used, currently available methods and apparatus for diagnosing 
AMI are the standard 12-lead electrocardiogram (ECG) and biochemical cardiac markers, such 
as creatine kinase isoenzyme-cardiac muscle sub-unit (CK-MB), tropinin T and I, and 
myoglobin. Less frequently used diagnostic techniques include computerized 12-lead ECG 
analysis, nuclear imaging (Technetium-99m-sestamibi), echocardiography, ECG exercise stress 
testing, continuous 12-lead ECG monitoring, and nonstandard ECG leads (i.e. right precordial 
chest leads). Despite this large number of option, standard 12-lead ECG and cardiac marker 
analysis remain the current mainstays of AMI diagnosis in the emergency room, likely due to 
availability and ease of use. 

[0007] The standard 12-lead ECG and cardiac markers each have several shortcomings in their 
abilities to quickly and accurately diagnose AMI. For example, 12-lead ECG analysis is a 
relatively fast diagnostic method but generally has low sensitivity and specificity in detecting 
AMI, particularly in patients with an acute posterior or lateral ML On the other hand, the most 
commonly used cardiac markers (i.e. CK-MB and the cardiac troponins) are not released until 
several hours after the AMI has occurred and, thus, do not provide for rapid diagnosis. Neither 
the 12-lead ECG nor the cardiac markers typically provide information regarding the location or 
size of an AMI. 

[0008] In an effort to enhance diagnosis of AMI, alternative ECG techniques have been 
proposed. Some techniques, for example, use greater than the standard number of 12 leads. For 
example, the PRIME ECG™ device, available from Meridian Medical Technologies, Inc. 
(Columbia, MD), includes 64 anterior leads and 12 posterior leads disposed in an electrode vest. 
Such ECG devices cover a larger area on the chest than conventional 12-lead ECGs, which 
typically improves their ability to accurately detect AMI. As with conventional ECGs, the 
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multiple-lead ECGs provide rapid diagnosis, as opposed to the slower diagnostic techniques 
involving cardiac markers. 

[0009] Although multiple-lead ECGs represent an improvement over conventional ECGs, they 
still have certain drawbacks. Currently available ECG systems, for example, typically detect 
AMI by comparing the ECG tracings of a patient to a preset baseline tracing. If the patient's 
tracing differs from the baseline by more than a specified threshold level, then the patient is said 
to be experiencing an AMI. Unfortunately, this method of diagnosis can at times lack sensitivity, 
specificity or both. This leads to some patients with AMI being discharged from the hospital 
without treatment, because their AMIs are not detected, while other patients are treated for AMI, 
for example with an angioplasty procedure, when they in fact did not have a AMI. Additionally, 
although multiple-lead ECGs are able to detect AMI accurately in many cases, they do not 
provide information regarding the location of the AMI on the heart, the size of the AMI, or 
similar diagnostic information which would help a cardiologist, surgeon or other treating 
physician to pinpoint and quantify the AMI. 

[0010] For the foregoing reasons, it would be advantageous to have methods and apparatus for 
enhancing diagnosis of AMI. More specifically, it would be advantageous to have improved 
methods and apparatus for detecting AMI, determining a location of an AMI on a heart, 
determining the size of an AMI, or any combination thereof. 

BRIEF SUMMARY OF THE INVENTION 
[0011] The present invention generally uses electrocardiogram data to enhance diagnosis of 
acute myocardial infarction (AMI). Typically, cardiac data from a patient, acquired from 
multiple-lead ECG devices (for example, between about 30 and about 130 leads), is compared to 
stored cardiac data from multiple AMI patients. The stored data also includes information such 
as presence of an AMI, the location of the AMI on the heart, the size of the AMI, and/or the like. 
The new patient's ECG can be matched with stored cardiac data that is similar to that of the new 
patient and, thus, the presence, location, size, and/or the like of the new patient's AMI can be 
determined. 

[0012] Generally, data is acquired from multiple patients, sorted, and stored, so that 
subsequent patients may be diagnosed based on the stored data. For example, cardiac data for 
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each patient may include signals from each of multiple ECG leads. The data from the leads are 
evaluated and reconstructed to create an integral map and a potential map by evaluating a 
segment, such as the ST segment, of all the signals and plotting them in a grid or similar matrix 
to create a unique cardiac map for each patient. The data for multiple patients is then sorted, so 
that patients with similar cardiac data maps are grouped together. The maps are also correlated 
with locations of ischemia in the patients' hearts, the locations being determined using standard 
methods such as thallium scans, angiography, and/or other diagnostic tests. When a new patient 
is evaluated via the multiple-lead ECG, the presence, location arid/or size of any AMI in the 
patient's heart may be diagnosed by comparing the patient's map to the stored data. 

[0013] In one aspect, a method for compiling data to enhance diagnosis of a myocardial 
infarction includes measuring cardiac data of multiple patients, using a multiple-lead 
electrocardiogram, sorting the cardiac data into groups, each group comprising cardiac data from 
one or more patients who experienced ischemia in a similar location on their hearts, and storing 
the sorted cardiac data. Optionally, the measuring may include acquiring electrocardiogram data 
from the patients with the multiple-lead electrocardiogram while the patients are suspected to be 
experiencing ischemia, confirming that the patients experienced ischemia, and identifying at 
least one location of the ischemia on the heart of each patient. Alternatively, the measuring may 
include acquiring electrocardiogram data from the patients with the multiple-lead 
electrocardiogram while the patients are experiencing temporary ischemia during percutaneous 
transluminal coronary angioplasty, and identifying at least one location of the temporary 
ischemia on the heart of each patient. In some embodiments, the identified location of ischemia 
on the feeaftehgM is taken from a group consisting of an anteroseptal location, an inferior 
location, a posterolateral location. 

[0014] In another aspect, a method for enhancing diagnosis of a myocardial infarction in a 
patient includes matching cardiac data acquired from the patient with stored cardiac data, 
wherein the stored cardiac data includes data from multiple ischemia patients sorted into groups 
of patients based on locations on the ischemia patients' hearts where ischemia occurred, and 
identifying at least one location of ischemia on the patient's heart, based on the stored cardiac 
data that matches the acquired cardiac data. In various embodiments, the method may also 
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include determining whether a myocardial infarction has occurred in a patient, determining a size 
of the myocardial infarction, and/or otherwise quantifying the myocardial infarction. 

[0015] In yet another aspect, a method for enhancing diagnosis of a myocardial infarction in a 
heart of a patient includes acquiring cardiac data from the heart using a multiple-lead 
electrocardiogram device, matching the cardiac data to stored cardiac data for at least one 
ischemia patient, wherein the location of ischemia on the ischemia patient's heart is known, and 
displaying information regarding at least one location of ischemia on the patient's heart, based 
on the matching stored cardiac data for the ischemia patient. 

[0016] A system for using stored electrocardiogram data to enhance diagnosis of a myocardial 
infarction in the heart of a patient may include a data storage module for storing 
electrocardiogram data of multiple patients, the patients having either experienced a myocardial 
infarction or temporary cardiac ischemia during percutaneous transluminal coronary angioplasty 
and the data being sorted according to locations of ischemia on the hearts of the patients, and 
computer software for enabling comparison of new electrocardiogram data to the stored 
electrocardiogram data to determine a location of ischemia in the heart of the patient. The 
system may optionally include a disposable panel including multiple electrocardiogram leads for 
acquiring cardiac data from a patient. The system may also optionally include a display module 
for displaying information related to a patient's heart. 

[0017] An apparatus for enhancing diagnosis of myocardial infarctions may include a database 
with stored electrocardiogram data from multiple patients with cardiac ischemia in known 
locations on the heart, the stored electrocardiogram data being sorted into multiple groups based 
on the known locations. Alternatively, an apparatus may include a disposable panel or similar 
device containing multiple electrocardiogram leads for measuring cardiac data in a patient. In 
other embodiments, an apparatus includes a multiple-lead electrocardiogram panel which may be 
coupled with a database for diagnosing and localizing ischemia and AMI in a patient. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[001 81 Fig. 1 A schematically illustrates a sensor system ha ving an array of sensing locations 
distributed across a patient's torso. 
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[00191 F ig. IB graphically illustrates the method for calculating an integral value across a 
selected time portion of a heart signal cycle from a single sensor location. 

100201 Fig. 1C illustrates a plot of a data matrix generated bv mapp in g the integral values with 
positions corresponding to the locations of the sensors across the patient's torso. 

100211 Fig. 2 schematically illustrates a method and computer program for separating a portion 
of an electrocardiogram signal relevant to a first portion of a heart from signals relevant to 
alternative portions of the heart. 

[00221 Figs. 3A-D graphically illustrate signal port ion se paration method steps employed bv 
the program and method of Fig. 2. 

[00231 Fig. 4 A graphically illustrates a database of known atrial paced heart cycles as 17 mean 
P wave integral maps. 

f00241 Fig. 4B and 4C illustrate 17 known right atrial ectopic origins associated with the 17 
mean P wave integral mans of Fig. 4A. 

r00251 Figs 5A-C illustrate a data base of mean P wave integral maps and associated locations 
of ectopic origins of the left atrium. 

10026] Figs. 6A-F illustrate correlations between integral maps of paced heart signal cycles 
obtained in different patients at a common region in the left atrium. 

[00271 Fig. 7 illustrates a database of ORS integral maps and associated ectopic origins within 
the right ventricle. 

10028] Fig. 8 illustrates a database of ORS integral mans and associated ectopic origins within 
the left ventricle. 

100291 Fig. 9 illustrates a diagnosis and treatment methodology using an arrhythmia location 
database and signal separator to help locate and treat the origins of focal AFib and other 
arrhythmias. 

100301 Figs. 9A and 9B illustrate arrhvthmo genie regions of the left and right ventricles. 
respectively. 
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[0031] Fig. 10 schematically illustrates a system and/or kit for diagnosing and/oi 
AFib and other arrhythmias, according to the principles of the present invention 



[00321 Figs. 1 1 A-l 1H illustrate exposed and skin-engaging surfaces of four panels supporting 
heart cycle sensors in an exemplary sensor array structure. 

[00331 Figs. 12 A and 12B are computer screen prints of an AFib localization program showing 
the localization output. 

DETAILED DESCRIPTION OF THE INVENTION 
fOOt^f -IflflML The present invention generally provides for enhanced diagnosis of myocardial 
infarctions, especially acute myocardial infarctions (AMIs). More particularly, the invention 
provides methods and apparatus for helping medical professionals or other users to detect an 
AMI, identify a location of the AMI on the heart, and/or determine the size of the AMI. From 
such determinations, a physician can make a more accurate assessment of which coronary artery 
is occluded and where treatment should be directed. 

fQ01jH -f00351 Most embodiments of the present invention employ multiple-lead 
electrocardiogram (ECG) devices and/or data from such devices to achieve such enhanced 
diagnosis. In some embodiments, for example, ECGs with between about 30 and about 130 
leads are used to compile a body surface map. One type of ECG device which may be used is 
described in PCT Application No. WO 01/67954 Al, (hereinafter "'954 application") previously 
incorporated by reference and attach e d h e r e to as App e ndix A . Such a device generally includes 
multiple panels, each panel including multiple ECG leads, which are configured to be placed on, 
under or around a patient's torso to collect cardiac data. (See, for example, Figs. 1 1 A-l 1H and 
accompanying description in '954 application.) 

[0036] The '954 application is directed to techn iques for localizing and/or treating atrial 
fibrillation and other arrhythmias. As explained in the f 954 application and illustrated in Fig. 1 A, 
the techniques of the present invention w ill gener ally make use of an array 10 of sensors 12 
distributed across anterior and posterior skin surfaces of torso T on patient P. Array 10 provides 
multilead electrocardiogram (ECG) data at a plurality of sensing locations distributed across 
torso T. typically at ov er 20 sensing locations, more preferably at over 40 sensing locations, and 
ideally at 62 or more sensing locations. 
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r00371 Sensors 12 generally comprise unipolar or bipolar electrodes coupled to the patient's 
skin, or to an alternative accessible body surface (Tor example via a transesophageal approach") 
suitable for measuring electrical body surface potential. Suitable electrode structures mav 
include those described in U.S. Patent Nos. 531 1.873 and 5.634.496. previously incorporated 
herein bv reference. Exemplary arrays for use in locations having large amounts of 
electromagnetic noise (such as an electrop hvsiolo gv lab or other location in which electrosurgery 
or electrical stimulation of tissues for intracardiac pacing is performed) was described bv 
Metting van Rim. A.C. et al. in IEEE Trans. Biomed. Eng.. BME-40:302-308: (1993V 
Alternative sensor array structures and associated data acquisition and manipulations 
components were described bv SippensGroenewegen. A. et al. in an article entitled. "Body 
Surface Mapping Purine Pacing at Multiple Sites in the Human Atrium: P wave Morphology of 
Ectopic Rizht Atrial Activation Circulation. 97:369-380 (1998): and bv Linnenbank. A.C. in a 
1996 thesis for the University of Amsterdam entitled. "On-Site Recording Analysis, and 
Presentation of Multi-channel ECG Data ". 

[00381 Referring now to Fig. IB. ECG data is preferably acquired simultaneously from each 
sensor 12 of array 10 at a sampling ra te of over about 500 Hz . ideally at a samp ling rate of about 
1.000 Hz or more. In so me embodiments, sequential sampling of sensor 12 from array 10 mav 
alternatively be used, and higher or lower sampling rates are als o feasible. When a lower 
sampling rate is used, the data mav be upsampled using multi-rate filter banks. 

f00391 Preferably, signals which are absent, for example, due to electrode obscurement bv 
defibrillator patches or lead dislodgment. may be deleted. Poor quality s ignals ma y also be 
visually and/or automatically identified and rejected. Such rejected signals mav be replaced 
using inte r polation of adjacent lead recording data. Interpolation techniques mav also be utilized 
to correct for offset variation among electrodes, and for linear baseline drifting. 

[00401 Graph 14 includes an ECG signal tracing 16 representing the variation in voltage over 
time, as sensed bv sensors 12. optionally at about 1 to 2 ms intervals. Signal tracin g 16 may be 
used to evaluate heart cycle signals from the heart of patient P. In general, one or more reference 
heart cycles will be selected for manipulation and comparison. The reference heart cvcle will 
typically be the heart cycle co incidin g with initiation of the arrhythmia for localizing focal AFib r 
or a p rem ature heart cvcle that has the same morphology as the cvcle that led to initiation of 
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focal AFib. ECG Tracing 16 can be used to determine a beginning 18 and end 20 of a time 
portion 22 of the reference heart signal cycle which is of particular interest for evaluating one or 
more regions of the heart. In the example illustrated in Fig, IB. a P wave onset may be 
determined bv the time at which the vo ltage progresses beyond 30 |uV while termination of the P 
wave may be defined at the atrial J-p oint. as is generally understood in the field of 
electrocardiography. Alternative criteria for P wave onset and offset might also be utilized, and 
automate d detection of time portion 22 is also feasible. Alternative time portions mav also be 
selected. 

[00411 Referring to Figs. IB and 1C. measurements made at each sensor 12 are preferably 
mapped onto a data matrix 24 according to the locations of the associated sensor. In the 
exemplary embodiment, a P wave integral numerical value 24 mav be calculated based on heart 
cycle signals 16 within selected time portion 22 for a particular sensor location N10. This 
calculated P wave integral va lue reflec ts the time/amplitude area of ECG signal at that sensor 
location within the selected time portion. Similar integral values are calculated for each sensor 
location, and the sensor values are mapped within data matrix 24 continuou sly from a portion of 
the data matrix associated with a front F of torso T. across a side of the p ati ent P. and to a back B 
portion of torso T. As shown in Figure 1C. the d ata matrix will often be presented graphically bv 
calculating lines of constant integral valu es 28 bas ed on the individual discrete integral values 
and their associated positions within the data matrix. In some embodiments, this information can 
be summarized bv presenting a single line 30 of zero integral value between a region of positive 
integral values 32 and a region of negative integral values 34. In much of the description which 
follows, the region of positive integral values 32 is presented as a shaded region within a 
graphically depiction of data matrix 24. Exemplary alternative data matrices may be presented 
with shad es of a fir st color (red, for example) for p ositive val ues, a second color (b lue, for 
examp le) for ne gative values, and optionally a third color (such as green) for zero. 

f00421 For localizing of certain arrhythmias, possibly including certain ventricular 
tachycardi as and some t ypes of atri al tachycardia, directly using measurements from sensors 12 
to calculat e integral values 26 for the selected time portion 22 mavbe sufficient to identify an 
arrhvthm ogenic region (which m ay be rel atively larg e) of a particular vent ricle, and in some 
cases, a particular atrium. Localizing directly from the sensed heart cycle signals is significantly 
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facilitated when the signals within the time portion of interest are predominantly indicative of 
activity within a candidate ectopic region of the heart. For example, when localizing ventricular 
tachycardia (VP. selecting a time portion dominated by the ORS complex in the signal can 
effectively localize arrhythmo genie foci or exit site, as more fully des cribed in the J. Am. Coll. 
Cardiol. 24:1708-1724 (1994). the full di sclosure o f which is incorporated herein by reference. 
This localizing of tachycardia foci within the ventricle mav be facilitated by the domination of 
the ORS complex in the signal of the overall body surface potential. 

[00431 Unfortunately, when localizing fibrillation foci within an atrium, the P wave (which can 
be indicative of activity within the at rium) will often be superimposed, either p artially or 
completely, bv the T-U wave. Physiologically speaking, the atrial activity of interest mav 
coincide with ventricular recovery of the preceding cardiac cycle. 

100441 To accurately localize focal triggers during the initiation of paroxysmal or persistent 
AFib. the present invention makes use of systems and methods for effectively separating a signal 
portion of interest from a superimposed signal portion, with the two signal portions often being 
separated from a single signal sensed from at least one single sensor location. These signal 
separation techniques are particularly advantageous when used to isolate the P wave from a 
simultaneously occurring T-U wave. It mav be possible in some circumstances to artificially 
separate these waves bv active overdrive pacing using an intracardiac catheter with a pacing 
period selected to avoid superimposition of these two signal portions during artificially initiated 
arrhythmia. For the reason described above, non-invasive electrocardiographic localization of 
atrial arrhythmias, particularly atrial premature b eats before or during invasive manning 
procedures is highly advantageous. As will be understood with reference to Figs. 2 - 3D. a 
ORST subtraction program helps to isolate and preserve the P wave morphology so as to enable 
trigger localization of focal AFib and other arrhythmias. The application of similar subtraction 
methodologies mav also enhance the ability of body surface mapping systems to localize 
triggers, exit sites, pathway insertion points, flutter waves, and/or fibrillation waves of other 
atrial arrhythmias such as atrial flutter, chronic AFib. and the like. ORST subtraction prog ram 
40 mav also enable application of inverse problem techniques to analysis of atrial arrhythmias, 
for example, when atrial depolarization is obscured bv the preceding ventricular repolarization. 
It should be understood that alternative signal separation methods and systems might also be 
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used, including those described in the following references, which are incorporated herein bv 
reference: Slocum. J. et al.. "Computer Detection of Atrioventricular Dissociation fro m Surface 
Electrocardiograms During Wide QRS Complex Tachycardias. " Circulation. 72:1028-1036 
(1985): Slocum. J. et al.. "Diagnosis of Atrial Fibrillation from Surface Electrocardiograms 
Based on Computer-Detected Atrial Activity. " J. ElectrocardiaL 25:1-8 (1992V Holm. M. et al. 
"Non-Invasive Assessment of the Atrial Cycle Length Purine Atrial Fibrillation in Man: 
Introducing. Validating and Illustrating a New ECG Method. " Cardiovasc. Res.. 38:69-81 
(1998): Bollman. A. et al.. "Frequency Analysis of Human Atrial Fibrillation Using the Surface 
Electrocardiogram and its Response to Ibutilide. " Am. J. Cardiol.. 81 : 1439-1 445 (19 98): and 
Ingemansson. M.P. et al.. "Autonomic Modulation of the Atrial Cycle Length bv the Head Up 
Tilt Test: Non-Invasive Evaluation in Patients with Chronic Atrial Fibrillation. " Heart. 80:71-76 
(1998). 

[00451 Referring now to Fig. 2. automated ORST subtraction program 40 uses an adaptive 
ORST template constructed from averaged ORST complexes combined with ECG body surface 
measurements to enable isolation of the otherwise obscured ectopic atrial activity. Generally. 
this approach allows the surface ECG measurem ents to retain their intricate spatial and temporal 
detail within the P wave morphology. Subtraction program 40 i s capable of unmasking and 
preserving subtle heart s ignal details of relatively low voltage P wave signal portions despite the 
obscuring superimposed relatively high voltage ORST complex. The ORST subtraction method 
of Fig . 2 i s described in more detail in U.S. Patent Application No. 60/189.513 filed March 15. 
2000. pre viously incorporated herein b v reference. 

[00461 As described above, the method of program 40 generally includes recording of unipolar 
ECG data from the array of torso sites in step 42. The measured signal will include both the P 
wave (which is of interest for AFib) and a superimposed QRST signal portion 44. 

[0047 ) In the exemplary embodiment, about 100 cardiac cycles of 62-channel ECG data are 
measured during sinus rhythm or atrial overdrive pacing. Op t ionally. ECG signals can be 
acquired during both sinus rhythm and atrial pacing. Fewer cycles may be used if the spatial and 
temporal variations of the ORST complex are relatively low. Typically, more than ten (10) 
cycles will be used, often more than 50 cycles for construction of the ORST template. As 
illustrated in Fig. 3 A. of f-line digital filtering o f the data designated for template construction 
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mav be performed using a 0,5 Hz high-pass filter, such as an IIR Chebvchev Tvpe-1 filter. This 
can hel p to correct for respiration-related baseline drifting. Additionally, a 100 Hz low-pass 
filter, such as a IIR Chebvchev Tvpe-1 filter, may be used to remove high-frequency signal 
artifacts. Additionally, a 50-60 Hz notch filter mav be used to remove line-frequency 
interference. Similar filtering mav be employed on the superimposed signals to be separated. 

[0048] Each filtered QRS complex 48 mav be identified using a complex-resonator/comb 
filter, together with a dual-edge threshold detection technique similar to that described bv Ruha. 
A. et al.. in an article entitled "A Real-Time Microprocessor QRS Detector System with a 1-ms 
Timing Accuracy for the Measurement of Ambulatory HRV'\ IEEE Trans. Biomed. Eng.. 
44:159-167 (1997V the disclosure of which is incorporated herein bv reference. Alternative QRS 
detection methods might also be used. 

[0049^ R wave fiducial points are marked and the average R-R interval is computed. The 
dominant QRS mor phology is identified, optionally using visual identification from pooled data, 
automated statistical methods, or the like. This dominant QRS morphology is used to select 
complexes for template creation. 

[00501 QRST template construction 50 mav be understood with reference to Fig. 3B. The 
selection of complexes for template creation mav be based on two criteria: QRS pattern, and R- 
R interval length. 

[0051 1 Regarding QRS pattern criteria, complexes of each cycle are compared with the 
template using parametric cross-correlations. In the exemplary embodiment, Pearson's 
coefficients are computed for a fixed time window sliding over a 20 ms time period. Once again, 
a variety of alternative cross-correlations mav be used. Each newly selected complex (for 
example, those having r > 0.98) is ali gned and averaged with previously selected cycles, and the 
QRS temp late is u pdated. Each complex which does not have adeq uate correlation with the 
tem plate is excluded, with these complexes often bein g ectopic or aberrantly conducted 
ventricular beats. 

[00521 Minimum R-R duration threshold is computed and only complexes having an R-R 
duration above the computed threshold are used to create the QRST template. For each complex, 
the O wave fiducial point is identified and the average Q-R interval is computed. The average 
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0-T interval (the length of the ORST complex) is computed optionally using a modified 
Bazett's formula. Additional correction for baseline drifting may be performed using linear 
interpolation after averaging a window prior to the O wave. An adaptive template may be 
constructed from the selected complexes by averaging their ORST intervals, optionally with the 
complexes aligned bv a window surrounding the R wave fiducial point. 

[00531 I n ste p 52. the above-described fiducial window in the ORS complexes, together with 
an additional fiducial time window around the peak of the T wave, are marked automatically. 
These fiducial windows are marked in both the ORST template, and in the superimposed signal 
48 containing a ORST complex together with a superimposed P wave. 

[00541 In step 54. the template 56 and the superimposed s ignals are ali gned. Typically. T 
wave window onset for the template is estimated as 0.64*OT. The T wave window is adjusted 
for each complex based on the OT interval length. Optionally, the operator can manually adjust 
any fiducial marker, such as by manipulating a mouse, joystick, bv putting a numerical value, or 
the like, the template's ORS fiducial points and T wave windows are aligned with the respective 
windows of the su perimposed ORST complex. A lignment of superimposed signal 48 and 
tem plate 56 mav be performed manually, or automatically by sliding windows over each other in 
10-ms increments and calculating cross-correlation coefficients. 

[00551 After alignment of the fiducial windows, the template is resampled as well as amplitude 
modulated 58. and then ORST subtraction is carried out as can be understood with reference to 
Fig. 3D and step 60. Template resampling and modulation are performed to compensate for 
discrepancies in duration (as a result of rate-related differences in the ORST int erval) and 
amplitude (as a consequence of variations in peak R-T wave voltage, which may be 
predominantly caused bv respiratory variation). This process is particularly geared toward 
obtaining optimal subtraction performance of the T-U wave complex. Optionally, cubic spline 
interpolation methods mav be used for template resam pling, although other multi-rate processing 
methods can also be used. 

[00561 The subtraction procedure from channel-to-channel mav be performed according to 
electrode position. This facilitates main tainin g sliding intervals as small as 10 ms while having 
the O-T interval dispersion accessed in multiple leads. Additional low-pass filtering can be 
carried out to smooth possible QRS leakage after ORST subtraction. The remaining ECG signal 
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after ORST s ubtraction features a P wave which is effectively is olated from the previously 
superimposed T-U wave. 

100571 After the ORST segment is effectively subtracted from each measured signal containing 
a P wave of interest, the morphology of the isolated P wave can be analyzed as described above. 
Specifically, an integral m ap may be computed of the separate d, previous ly superimposed P 
wave. This integra l ma p can be compared with a database of P wave maps created bv pacing, so 
that atrial tachycardia, atrial accessory pathway insertion sites, and focal triggers of paroxysmal 
or persistent AFib mav be localized using data from the surface ECG array 10. as shown in Fig. 
1 A. Advantageously, the method described herein above mav obtain high performance in the T- 
U wave range bv correcting specifically for differences in both the QRST duration, and in the 
voltage of the T wave. Additionally, the above-described QRST subtraction methodology makes 
use of a separate data set than the su perimposed wave to be separated, with the separate data set 
optionally being obtained during sinus rhythm or atrial overdrive pacin g to hel p insure that atrial 
and ventricular activity are clearly separated. As was also mentioned above, the superimposed 
wave or reference cvcle of interest will often comprise a single ectopic atrial heartbeat which can 
be readily se parated using the above-described method. 

[00581 Once the signal separation has been successfully completed for the reference heartbeat 
at each sensor location 12 of array 10. a P wave integral map mav be plotted from the data matrix 
26 as described above with reference to Figs. 1 A-C. Similar ORS integral data matrices and 
plots mav be generated for localizing ventricular arr hythmias. 

[00591 In many embodiments, data from such multiple-lead ECGs are used to create body 
surface maps, which are compared with stored body surface maps of multiple AMI patients. 
Again, for descriptions of similar body surface mapping methods and apparatus, in the context of 
cardiac arrhythmias, see the '954 application. In other embodiments, similar cardiac data is 
used, but the stored body surface maps are taken from patients who e xperiences experience 
temporary cardiac ischemia while undergoing angioplasty procedures. 

f0030 Hil060L Various embodiments of a method for enhancing diagnosis of AMI suitably 
include first collecting cardiac data. For example, in some embodiments, cardiac data may be 
acquired from patients undergoing percutaneous transluminal coronary angioplasty (PTCA). 
Cardiac data may be acquired using any suitable multiple-lead ECG device, but as discussed 
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briefly above, various embodiments employ one or more disposable panels, each of which 
contains multiple ECG leads. Disposable panels will generally be configured to be placed on, 
under or around a patient to measure ECG signals with the leads. When a PTCA patient is being 
treated, he or she will undergo a temporary cardiac ischemia when the angioplasty balloon or 
similar device is inflated. Cardiac data may be measured during that time, and a map made of 
that cardiac data may be correlated with the known area of ischemia. The area of ischemia is 
known, of course, because the surgeon or other physician conducting the PTCA procedure 
knows where the angioplasty balloon is located when it was inflated. 

f002j-HJQMlJ_In other embodiments, patients suspected of having an AMI could be assessed 
with an ECG device, such as the panels just described. The presence, location and/or size of the 
AMI could then later be confirmed by cardiac marker analysis, thallium scan, cardiac 
angiography, and/or other diagnostic methods. Angiography, echocardiography, and/or any 
other suitable studies may be used to provide further information about the AMI in each patient's 
heart. For example, the exact location of the AMI on the heart, the occluded coronary artery, the 
quantity of heart tissue affected by the AMI, and/or the like could be determined. This 
information could then be stored, along with ECG data, for each patient. In some embodiments, 
the information regarding the AMI and the ECG data will be processed to determine ECG 
characteristics that relate to AMI characteristics. For example, a particular ST elevation on an 
ECG may correspond to a particular location for an AMI in multiple patients. Such correlated 
information and data may then be used by medical professionals to more specifically diagnose 
AMI. 

f0023f IiMI62L In various embodiments, thenthe multilead ECG data will be used to create 
surface mapping of hearts for use in AMI diagnosis. Typically, such maps will be based on the 
following three sets of data or characteristics, though some embodiments may include less than 
all three: (1) infarct detection (rule in/rule out MI); (2) infarct localization (risk assessment to 
determine nature of acute treatment); and (3) infarct quantification (risk assessment to determine 
nature of acute treatment). Techniques for mapping are discussed in more detail in the "954 
application previously incorporated by reference. 

[00631 The '954 application is directed to techniques for localizing and/or treating atrial 
fibrillation and other arrhythm ias. As explained in the ! 954 application and illustrated in Figs. 
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4A-C. a graphical plot of a particular patient's P wave integral mav be used to localize an 
arrhvthmo genie region in an atrium bv comparing the P wave integral plot for the patient to a 
database of P wave integral plots having associated known ectopic foci within the right atrium. 
Each of the 17 plots of database 70 has an associated ectopic region (identified bv the encircled 
numbers illu strated in Figs. 4B and 4Q. the database having been gathered using pacing. 

[00641 The anterior-posterior view AP shown in Fig. 4B and the posterior-anterior PA view of 
Fig. 4C illustrate the right atrial cavity. Anatomical landmarks included in these diagrams 
include the superior vena cava SVC, and inferior vena cava IVC: the right atrial appendage 
RAA: the smooth right atrium SRA: the trabeculated ri ght atrium TRA: the crista terminalis CT: 
the fossa ovalis FO: the left atrium LA; the Eustachian valve EV: the coronary sinus os CSO: the 
tricuspid valve TV: the right pulmonary artery RPA: and the left pulmonary artery LPA. 
Methods for assemb ling a ri ght atrial database are described in detail in the J. Electrocardiol.. 31 
(Supp.^85-91 ( 1998V previously incorporated herein bv reference. A similar left atrial database 
80a. b. and c (generally referred to as left atrial database SO) is illustrated in Figs. 5A-C. Once 
again, the encircled numbers relate mean P wave int egral plots to specific endocardial regions of 
pacing shown in the anatomical diagrams. 

[00651 The mean P wave integral ma ps of left and right atrial databases 70. 80 feature extreme 
positions and zero line contours without positive and negative integral contour lines. Alternative 
plot formats, such as thr ee-dimensi onal or chest anatomv-based formats, map displays using 
various color schemes, and the like, mav also be used. 

[00661 The anatomical diagrams illustrated in Figs. 5 A-C present a posterior-anterior PA view 
and anterior-posterior AP view, and a left posterior-obliaue LPO view of the left atrium LA. 
Once again, major anatomical landmarks are highlighted including the left and right pulmonary 
arteries LPA. RPA: the superior and inferior vena cava SVC. IVC: the left atrial appendage 
LAA: the right atrium RA: the coronary sinus CS: the left ventricle LV: the left upper pulmonary 
vein LUPV: the left lower pulmonary vein LLPV: the right upper pulmonary vein RUPV: the 
right lower pulmonary vein RLPV: Bachmann's bundle BB: the mitral annulus MA: the anterior 
mitral valve leaflet AMVL: and the fossa ovalis FO. 

[00671 Referring now to Figs. 6A-F. six individual P wave integral mans included within 
group 6 were each obtained during pacing at the left upper and left lower pulmonary veins. 
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These six similarly located pacing sit es were grouped together within grou p 6 of left atrial 
database 8 0. and these plots w ere average d to produce the sixt h numbered mean plot of the left 
atrial database. The spatial compatibility of these patterns can be clearly seen, particularly with 
reference to the location and orientations of both the highest positive and negative integral 
values, as well as with reference to the zero line contour separating the shaded from unshaded 
regions. While each of these six patient-specific mans were generated using intracardiac pacing, 
naturally occurring ectopic origins may be identified bv comparing reference heart cycle signals 
measured during premature atrial beats, the onset of AFib. and/or atrial tachycardia (and 
o ptionally separated from superimposed signals described above) to these mean paced plots. 

K)0681 Referring now to Figs. 7 and 8 . a right ventricular database 82 and a left ventricular 
database 84 each include mean ORS integral mans for paced ectopic origins in the right and left 
ventricles, respectively. These ventricular databases are more fully described in an article by 
Peeters. H.A.P. etal entitled "Clinical Application of an Integrated 3-Phase Manpinz 
Technique for Localization of the Site of Origin of Idiopathic Ventricular Tachycardia 
Circulation 99:1300-131 1 f 1999) the disclosure of which is incorporated herein by reference. 

[0069] Referring no w to Fi g. 9. a localization/treatment method 100 may be performed bv first 
establishing a database of arrhvthmo genie regions and associated heart cycle signal 
characteristics, as described with reference to Figs. 4A-8 above. Body surface measurements are 
taken in step 102. typically using sensor array 10 d escribed with reference to Fig. 1A. As the P 
wave will often overlap other heart signal portions, localization of AFib and other arrhythmias 
will benefit from a signal separation pr ogram as described above regarding Figs. 2-3D. 

[00701 Once signal separation has been effected. P wave integrals for other selected indicia) 
are determined for each sensor 12. and data matrices assembled as described with reference to 
Figs. IB and 1C. Once again, it should be understood that these P wave integrals are calculated 
for a refer ence heart cycle signal of interest, usually for a premature atrial beat or for the heart 
cycle signal at the onset of an arrhythmia, for focal arrhythmias. Optionally, the data matrix mav 
be compared with the plots of the established database in step 104. This database comparison 
method m av simply involve visually selecting th e mean-pac ed integr al map which appears most 
similar to the plot for a specific patient. Alternatively, statistical correlation coefficients mav be 
generated bet ween the data matrix for the patient and each of the mean-paced groups of the 
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database. In some embodiments, different patient plots mav be prepared for comparison with 
different databases, for example, a P wave integra l ma y be calculated for comparison with atrial 
database groups, while a ORS integral plot mav be prepared for comparison with ventricular 
database groups. 

[00711 Once a mean-paced integral plot from the database has been selected as the closest 
correlation to the reference cycled plot for a particular patient, an arrhvthmo genie region 106 
associated with the corresponding mean-paced plot has effectively been identified. 
Arrhvthm o genie regions 106 associated with mean-paced plots 1-25 of left ventricular database 
84 and plots 1-13 of right ventricular database 82 are illustrated in Figs. 9A and 9B. respectively. 
In many embodiments, these arrhvthmo genie regions will be discrete locations based on the 
information within the as sociated da tabase. Preferably, arrhvthmo genie regions 106 will have 
surface areas of less than about 5 cm 2 . Optionally, the arrhvth mo genie re gions mav have an 
outer radius which is less than about 2.5 cm. ideally about 1.0 cm or less. In some embodiments, 
the arrhvthmo genie regions identified bv sensor array 10 on the patient's torso mav be small 
enough that no further localization is needed, and ablation of the ectopic site or e xit site within 
the arrhvthm ogenic regi on may proce ed without excessive col lateral damage. 

[00721 Referring now to Figs. 9 and 10, once an arrhvthmo genie region 106 has been 
identified, it will often be advantageous to further localize an ectopic origin or exit site 108 
within arrhvthmogenic region 106 using a pace m appin g catheter 110. Advantageously, this 
mapping of ectopic origin or exit site 108 may proceed rapidly within the limited confines of the 
arrhythmogenic region 106, thereby reducing fluoroscopy time and radiation exposure to patient 
and attending personnel, decreasing the trauma a ssociated with accessing alternative portions of 
heart H. and the like. In general, pace mapping is effected bv electrical stimulation of candidate 
ectopic origins within arrhvt hmo genic region 106 using a distal electrode pair of catheter 110. 
This can induce ectopic heart beats which can also be measured bv array 10 (or a similar array 
adapted for use in a high electromagnetic noise environment. Positioning of the catheter tip 
may be monitored using biplane x-ray imaging. 

[00731 The surface ECG corresponding to the paced beats is recorded, and the desired integrals 
and associated data matrix is generated, as described above. By comparing the data matrix plot 
of the induced arrhythmia to the database and/or original arrhythmia recordings, the focal origin 
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or exit site of the arrhythmia relative to the catheter position can be estimated, optionally using 
the meth od described in U.S. Patent No. 5.31 L873, previously incorporated here in bv reference. 
The catheter may be moved to the indicated alternative site, and the pacing and measurement 
steps repeated iterativelv until the ectopic site is found where the paced data matrix plot best 
correlates with the data matrix plot of the arrhythm ia. Once again, this iterative process is 
greatly ex pedited bv concentrating the ectopic origin or exit site search within arrhvthmogenic 
region 106 identified using sensor array 10. 

fQ02^ hJM74L Qnce maps are made for multiple patients, the maps may be sorted and stored. In 
many embodiments, for example, maps are generally stored according to the location, on each 
patient's heart, where ischemia occurred. In very general terms, when a new patient is then 
diagnosed using a multiple-lead ECG device, data from the new patient can be compared with 
data from previous patients to determine where the new patient's ischemia is located, whether 
the patient is actually suffering an AMI, the size of the AMI and/or the like. 

fQQ34 H00751 F or example, in some embodiments criteria for AMI detection will be developed 
based on the analysis of 62-lead ST integral or potential maps. Such ST maps will typically be 
developed using the presence of a specific amount of ST elevation (for example > 0.1 mV ST 
elevation) in a required minimum number of chest leads (for example at least 2 leads in the 
complete array). These data can be obtained by acquiring prospective surface ECG mapping 
data in a number of patients with AMI documented by cardiac marker analysis and/or coronary 
angiography. This will enable an assessment of the specific level of ST elevation and the precise 
minimum number of chest leads required to ensure the highest sensitivity and specificity for 
AMI detection. In addition, the information provided by the number of leads showing reciprocal 
ST depression during acute ischemia may also be utilized. As an alternative to the use of ST 
integral or potential maps, information provided in other segments of the ECG (i.e. QRS 
complex or QRST segment) may also be used. 

fQ025 H00761 hi many embodiments, a set of stored data such as database of 62-lead ST integral 
or potential maps will be developed based on a differentiation of the three major infarct 
regions/locations of the left ventricle which each relate to the major coronary arteries. These 
regions are: (1) anteroseptal MI (ASMI) due to an occlusion in the left anterior descending 
(LAD) coronary artery; (2) inferior MI (IMI) due to an occlusion in the right coronary artery 
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(RCA); and (3) posterior or lateral MI (PMI or LMI) due to an occlusion in the left circumflex 
coronary artery (LCx). The locations of these infarct regions can be verified using coronary 
angiography, for example. The classification in IMI may also include the presence or absence of 
right ventricular infarction (high-risk versus low-risk patient), which again may assessed using 
coronary angiography. Using methods of various embodiments of the present invention, infarct 
regions/locations will be specifically related to an area on the chest where a particular ECG 
change is measured. In some embodiments, for example, the infarct regions will be related to 
areas on a chest where ST elevation and/or reciprocal ST depression occurs during a specific 
time instant in a potential map or during a given time interval in an integral map. 

{QQ26 H00771 As described briefly above, in alternative embodiments, other methods are used to 
develop the above-described relationships between infarct regions and areas on the chest of ECG 
change. In some embodiments, for example, 62-lead ECG map patterns may be obtained during 
percutaneous transluminal coronary angioplasty (PTCA) performed in patients with coronary 
artery stenosis. In this approach, PTCA is be used as a surrogate model of AMI. At the time 
when a coronary stenosis is dilated during a PTCA, the vessel occlusion caused by balloon 
inflation results in temporary transmural ischemia which will generate an ECG pattern similar to 
the pattern obtained during ischemia caused by AMI. Acquisition of the ECG data during PTCA 
enables direct correlation of the map patterns with the angiographically validated balloon 
occlusion site. As an alternative to the use of ST integral or potential maps, information 
provided in other segments of the ECG (i.e. QRS complex or QRST segment) may also be used. 

f002?f IflflISL In some embodiments of the invention, information on the size of myocardial 
infarcts at each of the three major infarct locations (i.e. LAD, RCA, and LCx) will be obtained 
by developing criteria for analysis of 62-lead ST integral or potential maps that enable 
differentiation between patients with a proximal occlusion versus a distal occlusion of the 
infarct-related artery. In some embodiments, further differentiation may be made, for example 
by differentiating between proximal, middle, and distal occlusions. The location of a proximal 
versus a distal occlusion of the infarct-related artery may be verified using coronary 
angiography. As an alternative, PTCA may again be used as a surrogate model of AMI. 
Acquisition of the ECG data during PTCA would enable direct correlation of the map patterns 
with the angiographically validated proximal or distal balloon occlusion site. 
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fOO2^HQ022LIn various other embodiments, criteria may be developed for patients with non-ST 
elevation MI, patients with left (LBBB) or right (RBBB) bundle branch block, left ventricular 
hypertrophy (LVH), previous (i.e. old) myocardial infarction and/or the like. As with the criteria 
already discussed, these additional criteria may be determined by examining AMI patients via 
coronary angiography or by examining PTCA patients. Data from such examinations may then 
be used to derive criteria related to ECG data. 

fQQ2*H- [0080| Once data are collected and one or more sets of criteria for multiple-lead ECG 
diagnosis of AMI are developed, the data and criteria may be provided to medical professionals 
in any of a number of suitable forms, to enhance the diagnosis of AMI. In some embodiments, 
for example, a 62-lead ECG or other ECG device may be coupled with a database such that ECG 
data collected for a given patient may be quickly compared with data and/or criteria stored in the 
database. Such comparisons may enable the 62-lead ECG to detect an AMI, determine a 
location of an AMI, determine a size of an AMI and/or the like. 

[00811 Fies. 1 1 A-H illustrate skin-engaging and outer surfaces of four flexible senso r panels. 
These panel structures and the use of their array or sensors are more fully described in U.S. 
Patent Application No. 09/61 LI 79. filed on July 6. 2000. (now U.S. Patent No. 6.584343) the 
disclosure of which is incorporated bv reference in the '954 application, previously incorporated 
bv reference. 

f0030 Hflfl82L Typically, methods and apparatus of the present invention will provide for display 
of at least one diagnostic characteristic. For example, in many embodiments the detection of an 
AMI will be displayed, a location of the AMI will be displayed, and a size of the AMI will be 
displayed. In other embodiments, an occluded vessel may be listed or displayed and/or a section 
of the vessel which is occluded may be displayed. Display may be accomplished by any suitable 
means, such as an LCD screen, a computer screen, and the like. For example, display techniques 
such as those described in the -'954 application may be used. (See^ggg for example^ Figs. 12 A- 

r00831 T he '954 a pplication is directed to techniques for localizing and/or treating atrial 
fibrillation and other arrhythmias. As explained in the '954 application and illustrated in Figs. 
12A and 12B. screen prints 150a. 150b are at least a portion of the output of a localization 
computer program as described herein. Integral maps 152 are multi-color graphical 
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representations of th e mor phology derived from body surface potential maps 154 (also displayed 
in a multi-color format! Time portion 156 of reference cycle 158 is shown graphically, and may 
be separated from a superimposed signal bv inputting a "subtract" command 160, ideally using a 
gra phical user interface. Identified arrhvthmo genie sites 12 ar e out put with their associated 
probabilities, and are graphica ll y illustrated 164 relative to the adjacent anatomy. In the 
exemplary output, candidate sites of a selected database are also shown, with one or more of the 
identified sites being highlighted using color, blinking, an enhanced font or icon, or the like. 
Optionally, evele signals from different recorded segments 166 may be selected, and the cycles 
of interest mav be time scaled or zoomed 168 to show the desired evele intervals. Signals from 
selected le a ps or sensor panels may be displayed, and a variety of additional outputs may be 
provided. In the output illustrated in screen prints 150a. 150b. two different reference cycles 
from different arrhythmia events result in similar identified arrhythmo genie sites. Hence, the 
present invention can also provide confirmation of an arrhvthmo genie site identification from 
one or more confirmation refer ence cycle. 

f00841 In some embodiments, ECG readings of a patient include intuitive display formats in 
anatomical models of the heart. Existing imaging techniques may be used to verify the infarct 
location or infarct-related artery. Such techniques include, but are not limited to, coronary 
angiography (coronary artery anatomy assessment), nuclear imaging (regional perfusion 
assessment), and echocardiography (wall motion assessment). 

fOOa^ IMSSL Although the above is a complete description of a number of exemplary 
embodiments of the invention, other embodiments, as well as variations of the described 
embodiments, are also contemplated within the scope of the present invention. Therefore, the 
foregoing description is intended for exemplary purposed only and should not be interpreted tg 
limit the scope of the invention as it is described in the appended claims. 
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